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An implementation of a decoupled, single-input/single-output control approach for 
the Large Angle Magnetic Suspension Test Fixture is described. Numerical and experi- 
mental results are presented. The experimental system is a laboratory model large gap 
magnetic suspension system which provides five degree-of-freedom control of a cylindri- 
cal suspended element. The suspended element contains a core composed of permanent 
magnet material and is levitated above five electromagnets mounted in a planar array. 


1 INTRODUCTION 

A research effort is underway at NASA Langley Research Center to demonstrate the mag- 
netic suspension, positioning, and maneuvering of objects over wide ranges of attitudes. 
Future applications of this technology range from magnetic suspension of wind tunnel mod- 
els to advanced spacecraft experiment isolation and pointing systems. As part of this 
effort, a Large Angle Magnetic Suspension Test Fixture (LAMSTF) has been designed and 
built. The LAMSTF is a small scale laboratory model of a Large Gap Magnetic Suspension 
System (LGMSS) which provides five degree-of-freedom control of a cylindrical suspended 
element that contains a core composed of permanent magnet material. The suspended el- 
ement is levitated above five electromagnets mounted in a planar array. The LGMSS is 
a conceptual design of a ground based experiment which is to be used to investigate the 
technology issues associated with: magnetic suspension at large gaps, accurate suspended 
element control at large gaps, and accurate position sensing at large gaps [1]. The objec- 
tives of the LAMSTF effort were to investigate the feasibility of the LGMSS concept and to 
provide a test fixture for developing and demonstrating control approaches. A description 
of the LAMSTF and some of the control approaches which have been investigated are pre- 
sented in [2]. LAMSTF suspended element parameters and the field components generated 
by the electromagnets at the centroid of the suspended element are given in the Appendix 
to this paper. 

This paper presents the implementation of a single-input/single-output (SISO) control 
approach for the LAMSTF with numerical and experimental results. The control approach 
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which was implemented was developed in [3] using the extended linearized model developed 
in [4]. The control approach is proportional-derivative (PD), where the command torques 
and forces are functions of positions and derivatives of position. The design technique 
provides a dynamic compensator given the desired pole locations of the closed-loop system. 
Numerical results are obtained which compare the desired poles to the actual coupled pole 
locations for a candidate design. Experimental responses are also obtained on the LAMSTF 
testbed and compared to simulation results. 

2 EQUATIONS OF MOTION 

The equations of motion for the LAMSTF were developed in [4]. The LAMSTF configu- 
ration consists of five electromagnets mounted in a planar array. Figure 1 is a schematic 
representation of the LAMSTF configuration and defines the coordinate systems. The sus- 
pended element coordinate system consists of a set of orthogonaTx, y, z body-fixed axes 
which define the motion of the suspended element with respect to inertial space. The sus- 
pended element coordinate system is initially aligned with an orthogonal x, j/, z system fixed 
in inertial space. The open-loop equations of motion are determined by evaluating the forces 


z, z 



1 2 



Figure 1: LAMSTF Configuration 

and torques produced on the permanent magnet core by the magnetic fields [5]. They can 
be written as, 

X = f(X,I) (1) 


where 



X = [ 12 S 12* Vi Vy Vi 6y 

6 

and 


, n 


-2? 

-J? 

< 

II 

] 


In the state vector, 6 y and 6 Z are rotations about the y and z axis called pitch and yaw, 
respectively. The translations are x, y, and z, and 12 and V are time derivatives of the 


258 



corresponding positions. The input to the system are the five coil currents, denoted by 
I. In order to generate a linear model these equations are expanded about the nominal 
operating point X 0 ,I 0 using a Taylor series expansion and simplified using small angle 
approximations. Higher order terms are neglected in the expansion and motion about the 
uncontrolled axis (a:) is assumed to be zero. Details of the linearization are presented in 
[4]. The linearized equations have the form 


SX = ASX + BSI 


( 2 ) 


where 


A = 


dX 

dx 


and 


B = 


XoJo 


dX 

81 


XoJo 


Expanding these equations yields, 

ijly = vM x {-B x 9 y - B xz x- B yz y- B zz z)-vM x K z I 
Ictli = vM x (— B X 9 Z + B xy x + B yy y + B yz z) + vM x K. y I 
m c V x = vM x (-B xz 0 y + 2B xy 9 z + B xxx x + B xxy y + B xxz z) + vM x K xx I 

T1l C Fy VM X {By Z 6y + ( Byy B XX ) 9 Z + B X y X X + B X yy J/ + B X y Z Z ) + UMfK X yl 

™ C V Z = vM x (( B xx B zz ) 9y + B yz 9 z + B xzx x + B xzy y + B xzz z ) + vM x K. xz I 


( 3 ) 

( 4 ) 

( 5 ) 

( 6 ) 
( 7 ) 


where the B terms describe components and spatial derivatives of the magnetic field vector 
at the equilibrium point. The first subscript of B refers to a unit vector direction, while 
additional subscripts imply partial derivatives with respect to the coordinate system. The 
K coefficients are row vectors which define the fields produced by each coil per amp of 
current. I c and m c are inertia and mass of the core respectively, and vM x is the product 
of the core’s volume and magnetization. The first terms on the right in equations (3)-(7) 
are the torques and forces generated on the core due to perturbations in X , evaluated in 
the presence of the uncontrolled fields and gradients produced by the constant bias currents 
required to provide equilibrium suspension. The second terms are the torques and forces 
generated on the core by controlling the coil currents about the suspension currents. The 
controlled torques and forces can be written as, 


T c 

F c 


= BI 


( 8 ) 


where 


B = vMi 


-K z 

K y 

Kxx 

Kxy 

k tz 


( 9 ) 


For the LAMSTF configuration B is full rank 


and the currents required to produce given 
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command torques and forces become 


/ = B~ l 



( 10 ) 


For equilibrium suspension the torques and forces produced on the suspended element must 
be zero, except to counteract the effect of gravity. Therefore at equilibrium the suspension 


currents are 


I 0 = B 


-l 


0 
0 
0 
0 

L m c9 


( 11 ) 


Although B~ l decouples the five degrees of freedom in terms of force commands, the 
system dynamics are still highly coupled through the destabilizing bias terms in equations 
(3)-(7). The bias currents, I Q , are used to calculate the values of the bias fields and gradients 
presented in the Appendix. Many of these terms can be shown to be zero due to symmetry 
of the five-coil planar array. Referring to Figure 1, it can be seen that coil 1 is located 
symmetrically about the x,z plane and so cannot produce field components in this plane 
which are along the y axis. Furthermore, coils 2 and 5, and coils 3 and 4 form symmetric 
pairs about the i, z plane. These coil pairs have equal bias currents, therefore the projection 
along the y axis of fields from each pair is also zero throughout the x, 2 plane. This means 
that at X Q , the y component of the field is zero and all of its derivatives with respect 
to x and z are zero. Similar arguments can be used to show that B yy and B zz are also 
zero. Considering these zero terms, and making the additional approximation B xxz — 0 
and B xzx ~ 0, equations (3)-(7) become, 


I c tl y = —K B x 6 y — K B xz x + Ty C (12) 

I c tl- Z = -KB X 9 Z + T zc (13) 

m c V x = KB xxx x - KB xz e y + F xc (14) 

m c V y = KB xyy y + Fy C ( 13 ) 

m c V z = KB xzz z + F zc (16) 


where the K B terms are constants equal to the product of vM x and the corresponding field 
or gradient value evaluated at X 0 , / e . From these equations, it can be seen that the dynamics 
in y, z and 9 Z are uncoupled and can be analyzed as single degree-of-freedom systems. The 
dynamics in 9 y and x, however, remain coupled. The strength of this interaction and its 
effect on SISO control design will be addressed in the next section. 

The term KB X in equations (12) and (13) is negative and causes open-loop instability. 
These terms cause high frequency unstable modes referred to as compass needle modes. 
Compass needle modes occur because, with the LAMSTF configuration, in order to achieve 
gradients which generate the vertical suspension force, the core’s magnetization vector must 


260 



be aligned 180 degrees from the suspension field vector. These modes dominate the dynam- 
ics in pitch and yaw. The bias terms I\B XXX and KB xyy also cause unstable dynamics. 
These terms are similar to the unstable bias flux stiffness terms encountered with small gap 
magnetic bearings that use permanent magnet bias flux [6, 7]. The terms KB XZ and KB XZZ 
cause stable coupling between x and 9 y , and stable oscillations in z. Eigenvalues for the 
LAMSTF open-loop system are presented in Table 1. 


Mode Eigenvalue 


Compass needle 

-58.7793 

58.7793 


Compass needle 

-57.8061 

57.8061 


y Translation 

9.7764 

-9.7764 


Stable Coupling 

0.0000 

+ 7.9697i 


0.0000 

- 7.9697i 

2 Translation 

0.0000 

+ 0.9556i 


0.0000 

- 0.9556i 


Table 1: Eigenvalues of the open-loop system 


3 CONTROL SYSTEM EQUATIONS 

The control design technique allows the designer to directly place the poles of the closed 
loop system for each degree of freedom. Damping ratios and frequencies can be chosen to 
provide adequate response to disturbance inputs. The position of the suspended element 
is assumed to be known and is measured on the LAMSTF system by a set of five shadow 
sensors. As mentioned earlier, the control approach is PD, where the command torques 
and forces are functions of positions and derivatives of position. The command torques and 
forces can be written as 


Tyc ~ 

~{ P 6 y J r sRg 6 y 

(17) 

T- Zc = 

~ (Pe z + sRe z )0 z 

(18) 

F- - 

— ( P x + sR x )x 

(19) 

Fyc ~ 

~{Py + sRy)y 

(20) 

F ic = 

- (P z + sR z )z 

(21) 


Position and rate gains are denoted P and R, respectively, for each degree of freedom. 
Control of pitch rotation and x translation will be examined first since these are the only 
suspended element motions which are coupled. The approach is to close the loop around 
each axis independently and to determine the effect of the cross-coupling on the performance 
of the resulting system. Equations (17) and (19) can be written in matrix form as 
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(22) 


1 yc 

F xc 


e. 


= -Gf 

where Gf is the forward-loop transfer function matrix, 

Gf = 


Pe y + sRe y 


0 


0 


P x + sR x 


(23) 


Taking the Laplace transform of system equations (12) and (14) and putting them in matrix 
form results in 


' S 2 I C 9y ' 


s 2 m c x 



-I<B X -KB XZ 
— K B xz K B xxx 


9 


+ 


-t yc 

Fxc 


Substituting for T gc and F xc in (24) and collecting terms results in 

I c s 2 + R$ y s + P$ y + KB X A B xz 

I(B XZ m c s 2 + R x s + P x - K B xxx 



1 


* J 


The characteristic equation becomes 

/ c s 2 + Re y s + P$ v + A B x 


Vy v • * Vy 

KBr 


K Bxz 

m c s 2 + R x s + P x — KB X 


= 0 


Expanding the determinant yields 

( 


s 2 + ^y± + % + KB * 


2 RxS Px R Bxxx 
s H b — 


TTir 


m r 




( KB XZ ) 2 


m c 


(24) 


(25) 


(26) 


(27) 


where the system’s characteristic equation has been factored into two decoupled second 
order terms and a single coupling term. The coupling term in equation (27) is similar in 
form to the natural frequencies of the second order terms; however, it does not depend upon 
the feedback gains. It is possible to make this term negligible by increasing the position 
gains on the x and 6 y degrees of freedom. Ignoring the coupling term, the closed loop 
natural frequencies and damping for the 9 y mode can be written as 


= 


Pe y + KB X 


C&y 


Re., 


2 \/( P 0* + FB X ) J c 


(28) 


Similarly for the x control loop 


KB X 


W r = 


m. 


c« 


Rx 

2 v /(/ > r - K B xxx ) m c 


(29) 


Solving these equations for P and R yields design equations which allow for pole place- 
ment 


Pe u — F A B x 


Vy 

2 , 


P x — u 2 m c + I(B 


XXX 


Re y - 2 CdyUBylc 

R x — 2 Cx^x^c 


(30) 

(31) 
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The coupling term in equation (27) can only be ignored if the position gains are large, 
which implies that the SISO pole placement will be accurate only for sufficiently high fre- 
quencies. On the LAMSTF system this requirement was easily achieved. Table 2 shows 
design versus actual closed loop eigenvalues when the pitch and x loops are closed indepen- 
dently then analyzed as a coupled system. In each case the desired pole locations for the 
pitch and x modes were equal. Table 2 shows that the pole placement is inaccurate for low 
frequencies, but is reasonably accurate for frequencies above 75 rad/s, as expected. The 
table also shows very little variation between design and actual pole location for changes in 
damping. 


SISO design poles 
Frequency, rad/s Damping 

Coupled poles 

Frequency, rad/s Damping 

10.000 

0.707 

19.402 

22.117 

0.364 

1.000 

30.000 

0.707 

34.299 

24.971 

0.618 

0.849 

75.000 

0.707 

76.820 

73.133 

0.690 

0.725 

100.000 

0.707 

101.372 

98.608 

0.697 

0.717 

150.000 

0.707 1 

150.918 

149.075 

0.702 

0.711 

75.000 

0.100 

76.820 

73.133 

0.097 

0.102 

75.000 

0.300 

76.820 

73.133 

0.292 

0.307 

75.000 

0.500 

76.820 

73.133 

0.488 

0.512 

75.000 

0.700 

76.820 

73.133 

0.683 

0.717 

75.000 

0.900 

76.820 

73.133 

0.878 

0.923 


Table 2: Effect of coupling on accuracy of SISO pole placement 


The characteristic equations for the remaining degrees of freedom can be obtained in 
a manner similar to the pitch and x loops. Since these loops are uncoupled, the design 
equations are exact. The compensator parameters as a function of damping ratios and 
natural frequencies are given by, 


05 

1 

3 

It 

d? 

Re z 

— 2 Ce z u>e z Ic 

(32) 

P y = u>lm c + K B xyy 

Ry 

= 2 Q y u y m c 

(33) 

P z — oj x tti c + K B xzz 

R z 

= 2 ( z u> z m c 

(34) 


263 




4 IMPLEMENTATION 


The design method generates compensators which achieve approximate pole placement in 
the closed-loop system. The goal was to stabilize the experimental system and be able 
to maintain stability in the presence of disturbance forces. Poles were chosen to yield a 
system with a stiff response to disturbances and adequate damping to limit overshoot and 
oscillations. 

As a design example the natural frequency of each closed-loop pole pair was set to 
75 rad/s and the damping ratio of each was set to 0.707. Using equations (30)-(34) and 
parameters of the LAMSTF system, the position and rate gains were calculated. These 
gains are listed in Table 3. In the implementation of the PD controller it is desirable to 


Pe y 

= 4.94 

X 

10~ 2 

Nm/rad 

Re y 

= 5.84 x 10“ 4 

Nms/rad 

Pe, 

= 4.98 

X 

1(T 2 

Nm/rad 

Re z 

= 5.84 x 10" 4 

Nms/rad 

Pr 

= 1.25 

X 

10 2 

N/m 

R x 

= 2.35 

Ns/m 

Py 

= 1.25 

X 

10 2 

N/m 

Ry 

= 2.35 

Ns/m 

Pz 

= 1.25 

X 

10 2 

N/m 

Rz 

= 2.35 

Ns/m 


Table 3: Position and rate gains for example design 


limit the high frequency gain to minimize the effects of noise. Therefore, each PD loop was 
implemented as the following lead network, 


Pi 4 - sRj 
s/w r + 1 


(35) 


where the roll-off frequency, u T , is greater than the bandwidth of interest. In practice this 
was chosen to be 750 rad/second, an order of magnitude above the desired pole locations. 
A state-space model of the fully coupled system was developed and combined with the 
dynamic compensator to generate a continuous closed-loop system model. The damping 
ratios and natural frequencies of the poles of the closed-loop system are shown in Table 
4. The frequencies are higher than designed due to the effects of the roll-off pole in the 
implementation. 

The compensator designed above has been successfully implemented on the LAMSTF 
testbed. The suspended element position was derived from a set of five optical sensors 
arranged as shown in Figure 2. The sensors are based on power loss due to a shadowing 
of a collimated beam. The five measurements of the suspended element’s position were 
sufficient to calculate the position and orientation in five degrees-of-freedom. The sensors 
are accurate to about 10 microns and have a linear range of about ± 1 mm. Dynamics from 
the sensor electronics are negligible. 

An EISA-class 486 personal computer was used to implement the controller. Data acqui- 
sition, computation, and analog output were all handled by this computer. The controller 
was implemented as a set of discrete state space equations. The SISO continuous transfer 
functions were combined into a state space model and mapped into the discrete domain 
using a zero order hold transformation. The controller was implemented at a sample rate 
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Freq (rad/sec) 

Damping 

85.39 

0.648 

76.91 

0.730 

81.48 

0.703 

81.09 

0.669 

81.50 

0.703 

635.2 

1.000 

635.3 

1.000 

635.4 

1.000 

641.4 

1.000 

641.4 

1.000 


Table 4: Natural frequencies and damping of closed-loop poles 



Figure 2: Geometry of shadow sensor system 

of 1 kHz. A computational delay of 0.43 msec separated the input sampling and analog 
output times. Performance of the controller was demonstrated by subjecting the suspended 
element to equivalent pulse disturbance forces. Actuator currents corresponding to a given 
disturbance force or torque are calculated for each degree of freedom using B. These dis- 
turbances are implemented by adding these inputs to the closed-loop coil currents. Since 
the system has a nonzero steady-state error, the position of the suspended element tracks 
the input disturbance. Position of the element and the input torques and forces are plot- 
ted in Figure 3. The system remains stable and generally has a well damped response to 
the disturbances. Response in the yaw, x, y, and 2 degrees-of-freedom are all similar and 
consistent with design expectations. The response in pitch, however, contains slightly more 
overshoot and more oscillation than the other degrees-of-freedom. The cause of the under- 
damped response in pitch is currently being investigated and may be related to unmodeled 
dynamics from eddy current loops in the aluminum baseplate of the system. Additional 
testing and system identification are being performed to verify this hypothesis. A simula- 
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Figure 3: Experimental response to disturbance input 







tion of the closed-loop system was developed with MATLAB’s Simulink package [8]. The 
simulation models the continuous dynamics of the plant and considers the effect of sampling 
and computational delay in the implementation of the discrete controller. Figure 4 is a 
block diagram of the system as implemented in Simulink. 



Figure 4: Block diagram of Simulink system simulation 

Although the sampling frequency of 1 kHz is well above the closed loop dynamics, it 
was discovered that a pure analog analysis produced different responses than a simulation 
which considered the digital implementation. Figure 5 compares the full simulation, the 
experimental data, and an analog simulation. Response in pitch is underdamped experi- 
mentally, and the difference with respect to simulation can be easily seen. For the other 
degrees- of- freedom the digital simulation and actual response match quite well. The ana- 
log simulation, however, predicts a faster rise time and more overshoot in each case. A 
discrepancy between the frequency response of the continuous and discrete controllers was 
also noticed and appeared to be related to the high gain of the controller at the Nyquist 
frequency. Since the controller design was defined in the continuous domain, the differences 
between the digital and analog simulations are important to note. 


5 CONCLUDING REMARKS 

A decoupled control approach for a large gap magnetic suspension system has been pre- 
sented. The magnetic suspension system is a planar array of electromagnets which provides 
levitation and five degree-of-freedom control of a cylindrical permanent magnet. The con- 
trol approach assumes decoupled models for each degree of freedom. Position and rate gains 
for a dynamic compensator are computed based on desired pole locations. In the actual 
system, however, the system s dynamics remain coupled through bias terms resulting from 
the bias currents required to produce equilibrium suspension. The closed-loop performance, 
therefore, must be verified by applying the compensator to the coupled system model. 

This technique provides the control designer simple and intuitive parameters to adjust 
in order to achieve closed-loop performance. In order to investigate the effects of coupling 
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between pitch and x translation, an example design was performed using the parameters of 
the Large Angle Magnetic Suspension Test Fixture (LAMSTF). Results of the design indi- 
cate that the damping ratios and natural frequencies of the coupled axes differ only slightly 
from design values. The control approach has been experimentally demonstrated on the 

LAMSTF system. Transient responses to pulse inputs compare favorably with simulations 
of the closed-loop system. 
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A APPENDIX 


This appendix presents, in the form of tables, LAMSTF suspended dement parameters, 
electromagnet parameters, and components of fields and gradients (including second-order 
gradients) generated by the LAMSTF electromagnets at the centroid of the suspended ele- 
ment The LAMSTF contains a planar array of five room-temperature electromagnets with 
iron cores, mounted in a circular configuration. The configuration is shown schematically an 
Figure 1. For a more detailed description of the LAMSTF see [2]. The fields and grad, ents 
were calculated using VF/GFUN [9], including the pre- and post-processor OPERA with 
all iron cores modeled. Physical parameters of the LAMSTF are presented in table A . 
Electromagnet fields and first-order gradients generated by the suspension currents at the 
equilibrium point are presented in table A2. The fields, first-order gradients and second- 
order gradients generated by each coil at the equilibrium point are presented in table A3. 
It should be noted that only non-zero terms are included and the full set of components is 
not listed in the tables since Bij = Bji and Bijk = Bikj- 


Core diameter 

8.509 x 10 -3 m 

Core length 

5.08 x 10 -2 m 

Suspended element mass, m c 

22.124 x 10 -3 kg 

Suspended element inertia, I c 

5.508 x 10 6 

Core volume, v 

2.889 x 10 -6 m 3 

Core magnetization, M* 

7.785 x 10 5 A/m 

Suspension height 

0.1 m 

Electromagnet outer radius 

0.0825 m 

Electromagnet inner radius 

0.0475 m 

Electromagnet height 

0.105 m 

Iron core radius 

0.038 m 

Location radius* 

0.1375 m 

r * :i 


* Distance from center of array to axis of given coil 


Table Al: Physical parameters of LAMSTF 


Component 

Field strength, Tesla 

B x 

-8.1863e-03 

B xz 

9.6504e-02 

Bxxx 

4.9139e-01 

B X xz 

-2.4689e-04 

B X yy 

9.4051e-01 

B X zz 

-8.9865e-03 


Table A2: Values of bias fields and gradients at suspension point 
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Component 

Coil 1 

Fields, Tesla/Amp 
Coil 2 Coil 3 Coil 4 

Coil 5 

K x 

2.3100e-04 

7.1000e-05 

-1.8700e-04 

-1.8700e-04 

7.1000e-05 

Ky 

0 

2.2000e-04 

1.3600e-04 

-1.3600e-04 

-2.2000e-04 

K z 

-9.4000e-05 

-9.4000e-05 

-9.4000e-05 

-9.4000e-05 

-9.4000e-05 



First-order fie 

Id gradients, Tesla/m/ Amp 

Component 

Coil 1 

Coil 2 

Coil 3 

Coil 4 

Coil 5 

js 

XX 

2.1790e-03 

-1.5030e-03 

7.7200e-04 

7.7200e-04 

-1.5030e-03 

&xy 

0 

1.1960e-03 

-1.9360e-03 

1.9360e-03 

-1.1960e-03 

&xz 

-2.7230e-03 

-8.4100e-04 

2.2030e-03 

2.2030e-03 

-8.4100e-04 

Kyy 

-1.8920e-03 

1.7900e-03 

-4.8500e-04 

-4.8500e-04 

1.7900e-03 

Kyz 

0 

-2.5900e-03 

-1.6000e-03 

1.6000e-03 

2.5900e-03 

K zz 

-2.8700e-04 

-2.8700e-04 

-2.8700e-04 

-2.8700e-04 

-2.8700e-04 


Second-order field gradients, Tesla/m 2 / Amp 

Component 

Coil 1 

Coil 2 

Coil 3 

Coil 4 

Coil 5 

^ XXX 

3.4340e-03 

-1.8276e-02 

1.6559e-02 

1.6559e-02 

-1.8276e-02 

A xyx 

0 

-1.4560e-02 

1.3733e-02 

-1.3733e-02 

1.4560e-02 

TS 

f^XZT 

-5.3466e-02 

1.6371e-02 

-2.6790e-02 

-2.6790e-02 

1.6371e-02 

A X yy 

-2.1916e-02 

-1.1938e-02 

2.2896e-02 

2.2896e-02 

-1.1938e-02 

A xyz 

0 

8.7350e-03 

-1.4134e-02 

1.4134e-02 

-8.7350e-03 

TS 

xzz 

2.5400e-04 

7.8000e-05 

-2.0500e-04 

-2.0500e-04 

7.8000e-05 


Table A3: Field and gradient values for each coil, Tesla/ Amp 
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